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Abstract. A possible solution to the small scale problems of the cold dark matter (CDM) 
scenario is that the dark matter consists of two components, a cold and a warm one. We 
perform a set of high resolution simulations of the Milky Way halo varying the mass of the 
WDM particle (wwdm) an d the cosmic dark matter mass fraction in the WDM component 
(/w)- The scaling ansatz introduced in combined analysis of LHC and astroparticle searches 
postulates that the relative contribution of each dark matter component is the same locally 
as on average in the Universe (e.g. /w,0 = /w)- Here we find however, that the normalised 
local WDM fraction (/w,o / /w) depends strongly on to\vdm for niwDM < 1 keV. Using the 
scaling ansatz can therefore introduce significant errors into the interpretation of dark matter 
searches. To correct this issue a simple formula that fits the local dark matter densities of 
each component is provided. 
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1 Introduction 

Structure formation in the Universe is dominated by a mysterious dark matter component, 
whose energy density is about six times higher than the energy density made up of Standard 
Model particles. The kinds of structures which may form and their internal properties de- 
pends strongly on the as yet unknown physical nature of the dark matter particle. A common 
assumption is that dark matter is made of heavy, cold thermal relic particles that decoupled 
from normal matter in the very early Universe, generically referred to as cold dark matter 



Whilst there is a large body of indirect astrophysical evidence that supports CDM, 
especially on large scales, there are significant indications of possible shortcomings on small 
scales. First of all, CDM galaxy haloes contain a huge number of subhaloes [3-7], while 
observations indicate that only relatively few satellite galaxies exist around the Milky Way 
and M31 [8, 9]. One possible way out of this problem would be to have no or only very 
inefficient star formation in smaller CDM haloes [10, 11]. However, such models populate all 
the larger subhaloes with dwarf galaxies. Furthermore, the central densities of these largest 
CDM subhaloes seem to be about two times higher than the densities inferred form the 
observed stellar kinematics of the dwarf satellite galaxies [12-14]. 

Secondly, the inner density profiles inferred from galaxy rotation curves seem to be 
flatter than the NFW profile [15] used to fit simulated CDM halos (eg. [23, 24]). At the 
relevant scale of one percent of the virial radius the typical logarithmic density slopes are 
significantly steeper than the NFW form and there is substantial halo to halo scatter 7 = 
1.26 ± 0.17 (e.g. [18]). Fitting functions which are steeper than NFW at that scale provide 
better fits (e.g. [5-7, 17-19], which further increases the tension with the observations. 

Thirdly the observed number of dwarf galaxies in voids is smaller than expected from 
CDM [27]. 

One possible solution to these problems may be that dark matter is not cold, but 
warm (WDM). In this scenario dark matter particles are created relativistic unlike its CDM 



(CDM) [1, 2]. 
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counterpart and so 'free-stream' until relatively late times. Following [28], the free streaming 
length in the non-relativistic regime can be approximated by 



X FS ~ r H (t NR ) 1 + - log 



(1.1) 



tNR 



where tEQ denotes the epoch of matter-radiation equality after which the Jeans length starts 
to decrease. t^R is the time when the dark matter particles becomes non-relativistic. Due to 
the large mass of a cold thermal relic and its very early freeze out, the comoving horizon, r#, 
is extremely small compared with its warm dark matter counterpart, tNR strongly depends 
on the mass of the particle, that in turn therefore determines the characteristic mass scale 
below which primordial fluctuations are erased (e.g. [29]). 

There are several models motivated from particle physics that provide potential WDM 
candidates. The dynamics of the dark matter particle is mainly determined by its production 
history, i.e. if it once was in thermal equilibrium or not. The two most prominent candidates 
for these scenarios are the thermally produced gravitino [30] and the non-thermally produced 
sterile neutrino [31-33]. 

Refs. [34] and [35] recently reported a two sigma preference for an increased damping 
tail of the CMB power spectrum. As a possible explanation, they claim strong evidence for 
a fourth relativistic species. This is in line with measurements from atmospheric [36] and 
nuclear reactor [37] experiments, which find some evidence in favour of a sterile neutrino. 
Such light sterile neutrinos could be part of the cold + hot (rather than cold + warm) dark 
matter models. [38] have listed approximate conversion formulas for the mass of the warm 
particle for the most common production mechanisms, i.e. the thermal production of the 
gravitino [39], the resonant [40] and non-resonant production [31] of the sterile neutrino via 
oscillations with active neutrinos, and the sterile neutrino production via the decay of a scalar 



The power spectrum suppression of the latter two models is different from the usual thermally 
produced WDM and the above formulae give therefore only a rough estimate. These WDM 
candidates might help alleviate the small-scale issues of CDM, while maintaining the same 
large scale structure. However, the core created in such WDM models might not be large 
enough to match the observed inner density profiles of dwarf galaxies (e.g. [43, 45], and 
references therein). 

Another viable solution would be to have a mixture of CDM plus WDM - the C+WDM 
or mixed dark matter (MDM) model (e.g. [33, 46, 47]). Various theoretical motivations 
exist for C+WDM. The sterile neutrino for instance can be created via different production 
mechanism and appear in two disjunct, effective states, a cold and warm one (see [46] for 
details). One could also think of, besides the 'usual' cold neutralino candidate, to have a 
second dark matter particle, e.g. a thermal gravitino. 

Recently, a lot of effort has been dedicated to resolve the small-scale structure of CDM 
haloes [5-7] , however only a few simulations have been performed to study the WDM scenario 
[13, 38, 48], and a single one for the C+WDM model [49]. The aim of this work is to determine 



field [41, 42]: 
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the galactic dark matter halo properties in mixed dark matter cosmologies, especially the local 
dark matter densities of each component. A detailed analysis of substructure abundance and 
kinematics will be presented in a second paper (Anderhalden et al. 2012b, in prep.). 

This paper is organised as follows. In Section 2 the computation of the mixed dark 
matter transfer function is shown, as well as a summary of the numerical simulations. The 
internal structure of the Milky Way halo, i.e. density and phase space density profiles as well 
as the local dark matter abundance, is presented in Section 3. Implications of our results on 
dark matter detection are described in Section 4. 

2 Numerical Methods 

The transfer function of any non-standard cosmology can be written as the square root of 
the ratio of the power spectrum of cosmology i to that of CDM, i.e. 

The transfer function in a pure, thermally produced warm dark matter scenario is well fitted 
by [39, 50] 



T W dm(*0 = (l + (ak) 2u J (2.2) 
where v = 1.12 and a is a function of the warm particles mass, 

/ m x \ -i n / n x \ o-ii / h \ 1-22 . 

tt =°- 049 (i^) (oik) (tr) " Mpc - (2 - 3) 

The above form of the transfer function leads to a steep cutoff in the linear power spectrum, 
resulting in a strong suppression of small structures below the free-streaming mass (for a 
detailed discussion see e.g. [39, 46]). 

The qualitative behavior of the mixed C+WDM power spectrum differs from the pure 
warm dark matter case in one important aspect. Although the mass of the WDM particle 
still determines the cutoff scale, the presence of a cold component stabilises the drop-off and 
yields the transfer function to approach a constant plateau [46]. The size of this step is to 
first order completely described by the fraction of warm dark matter, i.e. 

f w — ^WDM _ ^WDM ^ 2 

^WDM + ^CDM ^DM 

Unlike Ref. [46] , we compute the transfer function of the cold and the warm component 
separately (even if they are not independent from one another, see Section 2.1). To illustrate 
this, a set of C+WDM linear power spectra at redshift z = 99 are plotted in Fig. 1. By 
keeping the mass of the warm particle constant (in this particular case mwDM = O.lkeV) 
and varying the fraction /w from zero (CDM) to one (WDM), one can see that the plateau, 
i.e. the step size, is dictated by the fraction, whereas the mass of the warm dark matter 
particle is responsible for the cutoff, i.e. the scale where the transfer function starts to 
differ from unity. Solid lines refer to the cold component of the mixed fluid, dashed lines 
to warm one. The higher the amount of warm dark matter, the more the plateaus of the 
two components approach each other, until finally a pure warm dark matter cosmology is 
reproduced (/w = !)• In the following we shortly describe how the aforementioned transfer 
functions are computed. 



-5/f 
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k [Mpc- 1 ] 

Figure 1. Evolution of the linear mixed dark matter power spectrum when keeping the mass of the 
warm particle fixed and varying the fraction /w- For reference the /w = (pure CDM) and /w = 1 
(pure WDM) cases are plotted too (black solid lines). In each case the colored solid lines refer to the 
cold component, whereas the colored dashed lines refer to the corresponding warm component of the 
mixed dark matter fluid. 

2.1 Mixed Dark Matter Transfer Function 

In order to run N-body simulations containing both cold and warm dark matter particles in a 
distinguishable way, the transfer function of each component has to be specified individually. 
On large scales, dark matter is very well described by a perfect fluid approximation. To study 
the behaviour of small perturbations embedded in a homogeneous and isotropic background, 
it is appropriate to linearize the complete set of equations (i.e. Continuity, Euler and Poisson 
equation (e.g. [51])). Since relativistic effects are negligible on scales below the horizon, we 
will restrict our computation to the Newtonian approximation at leading order, instead of 
solving the fully relativistic Boltzmann equation. 

The total dark matter density of the mixed (i.e. cold plus warm) fluid can be written 

as 

p(x) = p c (x) + pw(x), (2.5) 

where /? c (x) and /0 w (x) refer to densities of the cold and warm component respectively. Using 
the dimensionless density parameter, 5 = ^j=^, and defining fc = 1 — /w it follows 

Smdm = fw <5w + fc <>c- (2.6) 

By taking the ensamble-average over all /c-modes of 6, (\5k\ 2 ), the power spectrum of the 
mixed dark matter fluid can be entirely expressed in terms of /w / fc an d the CDM power 
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spectrum 1 , 



P 



MDM 



(*) = /w^WDM(fc) + JcPcDM(k) 

+ 2fwfcPcM k ) 

Jw r W + ic r C + 2 fwfcrwrc fbDM(fc) 



'MDM 



(k)P C Bu{k), 



(2.7) 



where we have introduced the ratios 



r c (*) - ^ aud rw(*) " M *° 



<5cDM 



^CDM 



(2- 



In order to set up initial conditions with a power spectrum for both the cold and the 
warm component individually, we also have to account for the interference term in Eq. (2.7). 
This can be achieved by solving the linearized Newtonian perturbation equations for a dark 
matter fluid consisting of a cold and a warm dark matter particle 2 (e.g. [51]): 



5c + 2-Sc 
a 

<5\y + 2 — <j\y 
a 



AttGpo 



fw <Hv + fc 8c 
fw $w + fc 5c 



(2- 



k 2 4 



5-w, 



where the derivatives are with respect to cosmic time and <to is the rms velocity dispersion of 
the warm component. Following [50], the velocity dispersion of the WDM particle at redshift 
z can be expressed as 



l + z 



/^WDM\I 

V 0.3 



h 
065 



keV 



4 
\ 3 



ra WDM ' 



km 



(2.10) 



and the rms velocity is do = 3.571i>o. gx in Eq. (2.10) is the number of degrees of freedom 
and can, for a generic thermal candidate, be set equal 1.5 [50]. Rewriting Eqs. (2.9) as a 
function of the scale factor, 



s'6 + Y J'c 



9vrGp 
9irGp 



fw <5w + fc 5c 
fw <5w + fc 5c 



(2.11) 



4 a? w ' 



a direct way to solve the above set of coupled differential equations is provided. The deriva- 
tives are now with respect to the scale factor. In the Newtonian approximation, dark matter 
perturbations do not have a significant growing solution before matter-radiation equality. We 
therefore choose a eq as the initial epoch for solving Eqs. (2.11) numerically. For adiabatic 
fluctuations, linear theory allows us to calculate the density contrast 5k(a eq ) to some time a 
(in our case the starting redshift of the simulation), via a transfer function T(k), 



5(k,aic) oc T(k)5(k,a eq ). 



(2.12) 



1 We again want to emphasise that this is only true in linear perturbation theory, neglecting all higher order 
interaction terms. 

2 For simplification we assume that the baryonic density perturbations follow those of cold dark matter. 
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Figure 2. Comparison between an empirical fitting function (solid black line, see Eq. (2.14)) for the 
plateau of the mixed dark matter transfer function and the one based on our two fluid approach (red 
squares) . 

It is straightforward now to calculate transfer functions for the cold and the warm component 
individually in such a way, that the power the of CDM model can be transformed into any 
mixed dark matter scenario: 



The above equation now provides an easy method to setup initial conditions for cosmological 
simulations. To test the validity of our simplistic approach, we make use of the fact that 
the mixed dark matter transfer function reaches a constant plateau at high values of k, fully 
determined by the fraction /w and compute T p i ateau at k = 10 3 h/Mpc for a set of different 
values of the warm dark matter fraction. Fig. 2 shows the comparison between our results 
(red squares) and an empirical fitting function of [46] (black solid line), describing a full 
Boltzmann code calculation for a mixed dark matter fluid. It is given by 



where g(ao) is a function of the energy density content of the universe (see [46] for details). It 
is evident that for small amounts of WDM the linearized approach slightly overestimates the 
exact value of T p i ateau , whereas it slightly underestimates Tpi a t eau when the warm component 
is dominant. However, this effect is extremely small, so it is sufficient to set up initial 
conditions for numerical simulations based on the procedure described above. 

Fig. 4 shows the effective initial power spectrum, calculated via Eq. (2.7), for all 
different cosmologies used in this work (see Table 1 for details). It is important to note 



P MDM( k ) = p CBM(k) ■ 2c iW 0) ^ -Pcdm(A;) ■ r£ jW . 



(2.13) 
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Figure 3. Density maps of the refined Milky Way haloes at redshift zero (R — r2oo)- From left to 
right: CDM, WDM 2keV and the most extreme mixed dark matter simulation f20 - O.lkeV. 

that Fig. 4 shows the power spectra for the mixed, i.e. cold plus warm fluid. Our initial 
conditions are set up by using the individual transfer function for each component separately 
(illustrated in Fig. 1). 

2.2 Simulations 

All numerical simulations have been carried out using the treecode PKDGRAV [52]. We per- 
form six different mixed dark matter cosmologies as well as a pure CDM and a pure WDM 
simulation with a particle mass of m w dm = 2keV. Our initial conditions are based on the 
WMAP7 cosmological model [53]: a 8 = 0.8, h = 0.7, fi dm = 0.227, O b = 0.046, A = 0.727, 
n s = 0.961 and are created with a parallel version of the GRAFIC package [54]. 

To begin with we run large scale simulations of two merged cosmological boxes with 
-^box = 40Mpc, all together using 2 x 256 3 particles. The merging of two boxes represents 
the mixture of the cold and warm component 3 and is accomplished by offsetting one of the 
boxes by half the mean interparticle separation in each direction, to ensure that particles 
do not lie on top of each other. The initial conditions for all eight simulations have been 
created using identical random seeds. We then select an isolated, Milky Way sized halo of 
mass M ~ 1O 12 M and re-run it at 8 3 higher mass resolution, leading to a dark matter 
particle mass of m p = 1.38 x 10 5 M and a physical gravitational softening of 355 pc in the 
refined region. The high resolution run was performed by a single level of spatial refinement, 
so in order to avoid mixing of heavy particles in the final halo, the Lagrangian region was 
chosen to be 3 x it^oo of the original selected object. The total number of high resolution 
particles within i?200 of the refined halo depends on the cosmological model, but is of order 
10 7 . Details of the simulations are listed in Table 1. 

The amount of cold and warm particles in the simulation is the same (2 x 256 3 ) and their 
masses are set by the fraction /c,w> i«e. they differ for /c,w / 0.5. In order to avoid mass 
segregation as much as possible, the two particle species have different softening lengths, 
proportional to the fraction /c,w [55], 

Cwarm = ^0 ' f^i C cold = ^0 ' fc^ i (2-15) 

3 We also choose that approach for the pure CDM & WDM simulations in order to have a one to one 
correspondence in the final halo. 
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Figure 4. Effective initial linear power spectra of all mixed dark matter fluids. 



Table 1. Details of the simulations. -R200 an d A^oo are measured with respect to 200 times the 
critical density. 



Label 


m WDM 


/w 


z=99 


-R200 


^200 




[keV] 


[%] 


[km s _1 ] 


[kpc] 


[xlO 6 ] 


CDM 










367 


9.88 


f05-0.05keV 


0.05 


5 


82.1 


362 


8.98 


f05-0.1keV 


0.1 


10 


32.6 


368 


9.46 


f20-0.1keV 


0.1 


20 


51.7 


365 


8.49 


f20-0.3keV 


0.3 


20 


11.9 


368 


9.20 


f50-0.3keV 


0.3 


50 


16.2 


367 


9.01 


f80-lkeV 


1.0 


80 


3.8 


368 


8.69 


WDM 


2.0 


100 


1.6 


367 


9.04 



where eo is chosen to be 1/50 of the mean inter-particle separation. Different gravitational 
softening lengths can in principle lead to a non conservation of momentum. However, PKD- 
GRAV solves this problem by a symmetrization of forces at all times. 

The internal structure of dark matter haloes is influenced by thermal motion of warm 
dark matter particles [45]. We add a thermal velocity component to all warm particles, 
computed as a function of the particle mass according to Eq. (2.10) [50]. We follow the 
procedure described in Ref. [45], i.e. creating a three dimensional set of Gaussian randomly 
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distributed velocities corresponding to the three velocity components of each particle and 
adding them to the Zeldovich velocities at initial redshift, z = 99. The variance of the 
distribution is dependent on the cosmology and is calculated via a(z) = 3.571 x v(z). 



3 The Inner Halo Structure 

An important characteristic of any non-standard cosmology is the possibility to obtain cored 
density profiles. In the following we present the measured densities as well as the pseudo 
phase space densities of the different Milky Way halo realizations. To illustrate the effect 
of cosmologies with non-zero free streaming, our refined Milky Way halo is shown in Fig. 3 
for three different cosmologies: CDM, WDM 2keV and our most extreme mixed dark matter 
model, i.e. f20 - 0.1 keV. The effect on the substructure number density is clearly visible. 

3.1 Density Profiles 

All density profiles of the simulations at redshift z = listed in Table 1 are shown in Fig. 
5. While the CDM halo shows the usual cuspy behavior in the central region 4 , a flattening 
becomes noticeable in all cosmologies containing a certain fraction of warm dark matter. The 
central densities in the innermost regions are damped by at least 20 per cent, even if the 
warm component is subdominant. In the f20 - 0.1 keV model a strong flattening is already 
evident at r ~ 10 kpc. 

In principle all mixed dark matter models can be brought into a hierarchical structure, 
ordered after its free streaming length. Considering a thermal relic with a Fermi-Dirac 
distribution and from the requirement to provide a correct dark matter abundance, it can be 
shown that the effective free streaming scale A is proportional to the fraction of warm dark 
matter times the average thermal velocity of the warm component [56]: 



The most extreme simulation according to Eq. (3.1), i.e. f20 - 0.1 keV, indeed shows the 
largest flattening of the inner density profile. However, the density profiles of all other 
cosmologies are not uniquely determined by Eq. (3.1), which should only be interpreted as a 
proxy for an effective free streaming of the combined, i.e. cold plus warm, mixed dark matter 



Fig. 6 shows the time evolution of the density profiles from redshift z = 1 until today. In 
all cosmologies, the halo has already fully virialized at redshift z = 1 and also the flattening 
of the inner density profile does not substantially evolve until z = 0. The profiles only show 
little evolution in the outer parts between redshifts z = 1 and z = 0, mostly due to infall of 
larger subhaloes. Fig. 7 shows the density profiles of each component. Fixing the amount 
of warm dark matter (/w)> an d at the same time varying its mass, creates large differences 
in the warm component. It differs by almost a factor of two at r ~ 1 kpc, while there are 
almost no differences the cold component. The outer part (r > 40 kpc) of the halo on the 
other hand is independent of the WDM mass and is mostly determined by the fraction /\\r. 

4 For completeness, we fitted the CDM run with three different fitting functions, of which the Einasto profile 
provides the highest accuracy. The best fit is skipped in Fig. 5 for the purpose of simplicity. 

5 This hierarchy is more apparent in the mass function, since the number density of collapsed objects is 
strongly dependent on the cutoff scale of the matter power spectrum, which in turn is a result of free streaming 
effects in the early universe. 
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Figure 5. The spherically averaged density profiles of the Milky Way sized halo in all eight cosmolo- 
gies at redshift z = 0. The data point represents the most recent value for the local dark matter 
density [57], i.e. p G = 0.008 ± 0.002 M pc" 3 . 



This is even true in the case for a warm dark matter dominated MDM model (see f80 - IkeV 
in Fig. 7). For comparison, Figs. 5 & 7 show a estimated local dark matter density [57]. The 
total dark matter density profiles of all models lie in the same range as the measurements of 
the local density. 

3.2 Pseudo Phase Space Density Profiles 

As pointed out by earlier works (e.g. [45, 58] and references therein), the formation of a core 
is related to the presence of a maximum in the phase space density distribution. It is useful 
to introduce the pseudo phase space density [59], Q = p/a 3 , where p is the dark matter 
density and a is the velocity dispersion. When Q is measured in CDM simulations, it usually 
follows a power Q(r) oc r~ a , with 1.8 < a < 1.95. Fig. 8 shows the pseudo phase space 
densitiy Q for all the mixed models as well as the pure CDM and the pure WDM case. Some 
of the models deviate from the power law behaviour (a ~ 1.9) only in the very central region 
of the halo. [45] showed that a core comparable with observations is only created in a pure 
warm dark matter scenario with a mass particle candidate, which is in conflict with recent 
Lyman-a forest data. Since we only ran a limited set of C+WDM simulations, we cannot 
tackle the question of whether or not such mixed dark matter models can be used to solve 
the cusp-core problem. 

Fig. 9 shows the time evolution of the pseudo phase space density component-by- 
component for our most extreme model, i.e. f20 - 0.1 keV. The black solid line represents Q of 
the cold component in the initial conditions (z = 99), the black dashed line the corresponding 
warm component. The black triangles on the other hand refer to the total pseudo phase space 
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Figure 6. Redshift evolution of the density profiles of all six mixed dark matter cosmologies. For 
a better understanding, only the density profile on top (f05 - 0.05kcV) is plotted in physical units 
whereas the subsequent profiles are each shifted by an additional factor of ten. 



density in the initial conditions. It is evident that the warm component shows little evolution, 
whereas the cold component drops by almost five orders of magnitude. 

3.3 Local Dark Matter Densities 

Unlike the C+WDM simulations in [46, 49], where every simulation particle represents the 
mixture of both species, our numerical approach allows us to investigate the composition of 
the local dark matter in a general C+WDM scenario. Let us define the local density ratio of 
either the cold or the warm component, 

Pc{r) + pw{r) 

Fig. 10 shows the contribution of WDM to the local dark matter density as a function of 
galactocentric distance, normalized to the universal fraction of warm dark matter, i.e. (/\y 
/ /w)- This quantity is found to be independent on /w, but to dependent strongly on the 
warm dark matter mass. The dashed vertical line in Fig. 10 further shows the position of 
the sun in the Milky Way and the black solid lines refer to an empirical fitting formula for 
the normalized local WDM fraction, given by 



nm ,r)=\l + ( .(^-) 2 { ] ^)}-\ (3.3) 
L VmwnM/ V r /J 



-mwDM' V r /. 

Using the Levenberg-Marquart method, we found that £ = 0.008 fits all C+WDM models 
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Figure 7. Density profiles of the mixed dark matter models, subdivided into the cold (solid) and the 
warm (dashed) components. For comparison, the most recent value for the local total dark matter 
density [57] is plotted. 

simultaneously with high accuracy. The deviations at small radii are probably due to higher 
order effects. 

We want to emphasize that Eq. (3.3) is explicitly valid only for a Milky Way size dark 
matter halo, since we expect the radial evolution of the warm/cold ratio to be dependent on 
the halo mass. In a pure warm dark matter cosmology, assuming isothermal density profiles, 
there is a dependence between halo mass and the expected core radius due to free streaming 
[44, 45], r c oc cr -1 / 2 oc M ha |^ 6 . We therefore expect the local density ratio at a fixed radius to 
be closer to the universal value when moving from Milky Way to cluster sized dark matter 
haloes. Additional simulations of dwarf (10 10 M ) and cluster (10 14 Mq) haloes are planned 
and will allow to quantify the halo mass dependence of these density ratios in detail. 

4 Implications for Dark Matter detection experiments 

Astroparticle searches for dark matter particles can be broadly divided into direct and indirect 
searches (see [60-64] for recent reviews). In this section, we discuss the implications of our 
results for both detection strategies. 

Direct dark matter searches aim to detect the recoil energy of detector nuclei struck by 
DM particles, assumed to be in the form of Weakly Interacting Massive Particles (WIMPs), 
a prototypical example of CDM candidates. The expected rate of events R is proportional 
to the product of the scattering cross section of WIMPs off baryons a cn times the flux of 
CDM particles through the detector, which in turn is proportional to the local dark matter 
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Figure 8. Pseudo phase space densities of all mixed dark matter cosmologies, as well as the pure 
CDM and WDM simulations. 

density, therefore 



In the literature, it is often assumed that pc,© = Pdm,©, he. that CDM completely saturates 
the measured DM density /9dm,© ~ 0.3 GeV cm -3 (see [21, 57, 65, 66, 68, 69] or a discussion 
of uncertainties on this quantity). This corresponds to the assumption that all the dark 
matter in the solar neighborhood is in the form of CDM particles, /c,© = 1- As we have 
seen, however, dark matter could in principle be made of a mixture of particles, in which 
case this assumption is explicitly violated. 

This is particularly important at a time when the Large Hadron Collider (LHC) is 
probing the parameter space of some of the most popular candidates for dark matter, such 
as the supersymmetric neutralino. In case of discovery, the question of whether the newly 
discovered particles compose all of the DM in the Universe can probably be answered only 
through a combined analysis that also includes direct and/or indirect detection experiments 



In order to combine LHC and direct detection data one can follow two strategies [70]. 
The first is to identify CDM with DM locally and on average in the Universe. Since the aim 
of the analysis is to prove that the new particles compose all of the DM in the Universe, this 
analysis is clearly not satisfactory and it can at best provide a consistency check for a DM 
interpretation of LHC data. 

To go beyond these simplistic assumptions, one can make the much less stringent ansatz 
fc,& = fc [70], which corresponds to assume that the fraction of CDM is the same locally as 



R OC (TcnPCe- 



(4.1) 



[63]. 
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Figure 9. Redshift evolution of the phase space densities for each the cold and the warm component. 
Shown is the example of our most extreme model, i.e. f20 - 0.1 keV Mixed Dark Matter. 



on average in the Universe. The rate of events can therefore be written as 

R OC (JcnPCQ = CTcnPefcQ = <?cnPQ)fc ( 4 - 2 ) 

Since this assumption proves to be very effective in breaking the degeneracies in the param- 
eter estimations of DM properties for combined LHC/direct detection searches, it is worth 
understanding the limits of its applicability to the problem at hand in light of the results 
discussed in the previous section. 

In the limit where the two (or more) components of dark matter are perfectly cold, one 
does not expect any deviations from the scaling ansatz above, since the two components are 
perfectly mixed in galactic halos down to the smallest scales. The results of our simulations 
allow us to precisely identify the regime where the fc( r )/fc starts to deviate from unity, i.e. 
when the scaling ansatz fails. Fig. 10 accurately provides this information, since it shows 
the behavior of this ratio as a function of the mass of the warm component and the distance 
from the galactic center. 

In Fig. 10 we observe deviations from perfect scaling at the solar radius by up to 30% 
for our models, although even larger deviations are possible (see Eq. (3.3)). Therefore, using 
the scaling ansatz introduces a significant systematic uncertainty and one should use Eq. 
(3.3) for a better estimate of the local densities. We also stress that the presence of a WDM 
component always leads to the inequality /c > /c- 

As for indirect detection, which is based on the observation of the annihilation or decay 
products of DM particles, we find as expected opposite effects on strategies for WDM and 
CDM searches. As we have seen, in fact, the joint evolution of these two components tends to 
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Figure 10. Local warm dark matter contribution normalized to the universal WDM fraction. The 
solid black lines show the fitting function of Eq. (3.3). The dashed line indicates the position of the 
sun in the Milky Way (r Q w 8 kpc). 

increase the relative contribution of CDM over WDM inside virialized structures. The specific 
enhancement (reduction) of gamma-ray (X-ray) fluxes expected from the annihilation (decay) 
of CDM (WDM) particles depends on the specific observational target and on the properties 
of the WDM component, and it is outside the aim of this paper. We limit ourselves here to 
notice that since all annihilation fluxes scale with the square of the number density of CDM 
particles, the predicted annihilation radiation depends strongly on the specific behaviour of 
fc( r ) in the virtualized halo under consideration. 

5 Summary and Conclusion 

We have performed a set of high resolution simulations of a Milky Way sized halo in different 
C+WDM cosmologies, varying the fraction and the mass of the warm component. For each 
run we therefore set up initial conditions, consisting of two distinct sets of dark matter 
particles, each species initialised with the corresponding modified matter power spectrum. 
This procedure enables us to track the density profiles of each component and to test the 
scaling ansatz, i.e. how the density ratio of the cold and warm components changes as a 
function of distance to the galactic center. We find: 

• The transfer functions of the cold and the warm component in a general C+WDM 
cosmology are strongly correlated. Unlike in pure WDM, where T{k) rapidly drops 
to zero after the cutoff, the presence of a cold component also stabilises the warm 
component and it approaches a constant plateau at high k. This leads to non negligible 
power in the warm dark matter component even at the smallest scales. 
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• Due to thermal velocities of the WDM particles, all mixed dark matter scenarios in 
this work show significantly reduced total dark matter densities in the halo center. The 
density profile of our warmest model, i.e. the f20 - 0.1 keV, shows a reduced density 
already at large radii (8 kpc) and a possible hint of a flat core in the very inner region 
(~ 1 kpc). 

• The assumption that the fraction of any component in a mixed dark matter model is 
the same locally as on average in the Universe, is violated in some of the scenarios. The 
presence of a warm dark matter component leads to the inequality /w < fw (/c — fc)- 



• The normalised WDM fraction (fw( r ) / fw) 1S proportional to the squared mass of 
the WDM particle and it is practically independent on /w- This is surprising, since 
mwDM is responsible for the cutoff of the linear C+WDM power spectrum, whereas 
the fraction determines the plateau at small scales. 

• In absence of further information on the nature of any dark matter component, our 
results on the local density introduce systematic uncertainties in the data analysis of a 
combined, i.e. accelerator plus direct detection, dark matter search. 

• For the interpretation of dark matter search data in generic mixed dark matter cos- 
mologies one should use Eq. (3.3) instead of the scaling ansatz, i.e. /c/w,© = fc/w- 
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